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ABSTRACT 

The exoplanets known as hot Jupiters — Jupiter-sized planets with periods less than 10 days — 
likely are relics of dynamical processes that shape all planetary system architectures. Socrates et 
al. (2012) argued that high eccentricity migration (HEM) mechanisms proposed for situating these 
close-in planets should produce an observable population of highly eccentric proto-hot Jupiters that 
have not yet tidally circularized. HEM should also create failed-hot Jupiters, with periapses just 
beyond the influence of fast circularization. Using the technique we previously presented for measuring 
eccentricities from photometry (the "photoeccentric effect"), we are distilling a collection of eccentric 
proto- and failed-hot Jupiters from the Kepler Objects of Interest (KOI). Here we present the first, 
KOI-1474.01, which has a long orbital period (69.7340 days) and a large eccentricity e = 0.81 jlo o?; 
skirting the proto-hot Jupiter boundary. Combining Kepler photometry, ground-based spectroscopy, 
and stellar evolution models, we characterize host KOI-1474 as a rapidly-rotating F-star. Statistical 
arguments reveal that the transiting candidate has a low false-positive probability of 3.1%. KOI- 
1474.01 also exhibits transit timing variations of order an hour. We explore characteristics of the 
third-body perturber, which is possibly the "smoking-gun" cause of KOI-1474.01 's large eccentricity. 
Using the host-star's rotation period, radius, and projected rotational velocity, we find KOI-1474. 01's 
orbit is marginally consistent with aligned with the stellar spin axis, although a reanalysis is warranted 
with future additional data. Finally, we discuss how the number and existence of proto-hot Jupiters 
will not only demonstrate that hot Jupiters migrate via HEM, but also shed light on the typical 
timcscale for the mechanism. 
Subject headings: planetary systems 



1. INTRODUCTION 

The start of the exoplanet era brought with it the dis- 
covery of an exotic new class of planets: Jupiter-sized 
bodies with short-period orbits (P < 10 days), com- 
monly known as hot Jupiters (Mayor fe Quelozl 119951 : 
iMarcy et al.l 119971 ). Most theories require formation of 
Jupiter-sized planets at or beyond the so-called "snow 
line," located at roughly a few AU0 and debate the 
mechanisms through which they "migrated" inward to 
achieve such small semimajor axes. The leading theo- 
ries fall into two categories: smo oth migration through 
the proto planet a ry disk (e.g. IGol dreich & Trcmainc 
19801 IWardl Il997l: lAliber t et al.l 120051: llda fc Linl Effll 
Bromley fc Kenvonll2011[ ). or what TSocrates et aTT(|2012l ) 



(hereafter S12) term high eccentricity migration (HEM), 
in which the planet is perturbed by another body onto 
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an inclined and eccentric orbit that su bsequently cir- 
cular i zes through tidal d i ssipation (e.g |Wu fc Murrayl 
Rasiol [2001 IFabryckv fc Tremaind 120071: 



20031: iFord fc Rasid 120061: IFabryckv fc 
Naoz et al.ll201lHWu fc LithwickboilT 



From the present-day orbits of exoplanets we can po- 
tentially distinguish between mechanisms proposed to 
shape the architectures of planetary systems during the 
early period of dynamica l upheaval. In this spirit, 
iMorton fc Johnson! (|2011al ) used the distribution of stel- 
lar obliquities to estimate the fraction of hot Jupiters on 
misaligned orbits and to disting uish between two specific 
migration mechan i sms (see als olFabrvckv fc Winn|[20ul: 
ITriaud et al1[20Tot IWinn et al.ll2010D : ? recently applied 
a similar technique to estimate the relative contributions 
of two different mechanisms. However, deducing dynami- 
cal histories from the eccentricity distribution of exoplan- 
ets poses a challenge because most hot Jupiters have al- 
ready undergone tidal circularization and "cold" Jupiters 
at larger orbital distances may have formed in situ. Fur- 
thermore, type-II (gap-opening) migration may cither 
excite or damp a planet's eccentricity through resonanc e 
torques (IGoldreich fc Sarilp o M ISari fc Coldreiclll27iol . 
Finally, iGuillochon et al.l ( 20111 ) find evidence that some 
hot Jupiters may have undergone disk migration either 
prior to or following scattering. In the latter case, disk 
migration may have damped their eccentricities. The 
eccentricity distribution is potentially shaped by a com- 
bination of HEM, tidal circularization, and planet-disk 
interactions. 

Motivated by the HEM mechanisms proposed by 
IWu fc Murray! (|2003| ) and others, S12 proposed an ob- 
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servational test for HEM. As an alternative to modeling 
the distribution of eccentricities, they suggested that we 
look for the individual highly eccentric, long-period pro- 
genitors of hot Jupiters caught of the act of tidal circu- 
larization. S12 identified HD 80606 b as one such pro- 
genitor, which was originally discovered by radial veloc- 
ity (RV) measur ements of its host star's reflex motion 
(jNaef et al.ll200lD and later found to transit along an or- 
bit t hat is misaligned with respect to its ho st star's spin 
axis (|Moutou et all 120091 : IWinn et all 12009( 1. From sta- 
tistical arguments S12 predicted that if HEM produces 
the majority of hot Jupiters, the Kepler Mission should 
detect several "super-eccentric" Jupiters with orbital pe- 
riods less than 93 days and eccentricities in excess of 0.9. 
A couple of these planets should be proto-hot Jupiters, 
with post-circularization semimajor axes in the region 
where all hot Jupiters have circularized (i.e. P < 5 
days). Several more eccentric planets should have final 
periods above 5 days, in the region where not all hot 
Jupiters have circularized; these planets may be "failed- 
hot Jupiters" that will never circularize over their host 
stars' lifetimes. A failed-hot Jupiter may have either 
halted at its post-HEM location due to the tidal circu- 
larization timescale exceeding the age of the system, or 
undergone some tidal circularization but subsequently 
stalled after a perturber in the system raised its peri- 
apse. S12's prediction is supported by the existence of 
super-eccentric eclipsing binaries in the Kepler sample, 
which are al so thought to hav e been created by HEM 
mechanisms (iDong et al.ll2012T l. 

To test the HEM hypothesis we are "distilling" ec- 
centric, Jupiter-sized planets from the sample of an- 
nounced Kepler can didates using the publicly released 
Kepl er light curves (jBorucki et al.l 120111 : iBatalha et al.l 
12012( 1. We described the distillation process and 
our technique for measuring eccentricities from tran- 
sit light curves based on th e "photoeccentric effect" 
in iDawson fc Johnsonl ()2012f) . hereafter Paper I. In 
summary, eccentric Jupiters are readily identi fied by 
their s hort ingress/egress/ total transit du rations ([B arnes 
20071: iFord fc Rasid 120081 iBurkd 120081 : iPlavchan etafl 
20121 Kane et al.ll2012h . A Markov-Chain Monte Carlo 
(MCMC) exploration of the posterior distributions of 
the transit parameters, together with a loose prior im- 
posed on the stellar density, naturally accounts for the 
eccentricity-dependent transit probability and marginal- 
izes over the periapse angle, yielding a tight measurement 
of a large orbital eccentricity (Paper I). 

Here we present the first eccentric, Jupiter-sized can- 
didate from the Kepler sample: Kepler Object of In- 
terest (KOI) number 1474.01. We find that this eccen- 
tric candidate also has large transit-timing variations 
(TTVs). In fact, the TTVs are so large that they 
were likely missed by the automatic TTV-detection algo- 
rithms, as they w ere not listed in a recent cataloging o f 
TTV candida t es (IFord et al.l [20121: ISteffen et al] 12012) . 
iBallard et al.l (|2011f ) recently deduced the presence and 
planetary nature of the non-transiting Kepler-19c from 
the TTVs it caused in the transiting planet Kepler-19b, 
demonstrating the viability of detecting non-transiting 
planet s through TTVs. More recently. iNesvornv et al.l 
(2012) characterized a Saturn-mass non-transiting planet 
using this technique. Thus the TTVs of 1474.01 may 
place constraints on the nature of an additional, unseen 



companion, thereby elucidating the dynamical history of 
the system. 

In Sj2l we present the light curve of KOI-1474.01. In 
Sj3j we characterize the host-star KOI-1474 using Ke- 
pler photometry, ground-based spectroscopy, and stel- 
lar evolution models. In <j4j we estimate the candidate's 
false positive probability (FPP) to be 3.1%. In we 
measure KOI-1474. 01's large eccentricity, investigate its 
TTVs and the perturbing third body that causes them, 
and measure the projected alignment of the transiting 
planet's orbit with the host star's spin axis. In J)J we 
place KOI-1474.01 in the context of known hot Jupiters, 
proto-hot Jupiters, and failed-hot Jupiters, and explore 
whether KOI-1474.01 is a failed-hot Jupiter that will re- 
tain its current orbit or a proto-hot Jupiter that will 
eventually circularize at a distance close to the host star. 
We conclude in ^7]by discussing the implications for plan- 
etary system formation models and suggesting directions 
for future follow up of highly eccentric planets in the Ke- 
pler sample. 

2. KOI-1474.01: AN INTERESTING OBJECT OF INTEREST 

KOI-1474.01 was identified bv lBorucki et all (|2011( 1 as 
an 11.3 Rq candidate that transit s its 1.23 M w , 6498 K 
host star every 69.74538 days (jBatalha et al] I2010I 1. 
With a Kepler bandpass magnitude Kp = 13.005, the 
star is one of the brighter candidates in the Kepler sam- 
ple, making it amenable to follow up by Dopplcr spec- 
troscopy. We retrieved the Q0-Q6 data from the Multi- 
mission Archive at the Space Telescope Science Institute 
MAST) and detre nded the light curve using AutoKep 
Gazak et al.l l2012). We identified eight transits (Figure 
[1]), which together reveal three notable properties: 

1. When folded at a constant period, the transits are 
not coincident in phase. Indeed, some fall early or 
late by a noticeable fraction of a transit duration. 

2. The transit durations are short for a planet with 
such a long orbital period (the total transit dura- 
tion, from first to last contact, is Tu—2.92 hours, or 
0.17% of the 69.74538 day orbital period). Yet in- 
stead of the V shape characteristic of a large impact 
parameter, the transit light curves feature short in- 
gresses and egresses — corresponding to a planet 
moving at 3 times the circular Keplerian velocity 
[based on the Kepler Input Catalog (KIC) stellar 
parameters] — and a nearly flat bottom, implying 
that either the planet has a large eccentricity or or- 
bits a very dense star (see Paper I) . T he candidate's 
repor ted a/R* = 129.0525 ± 0.0014 (jBorucki et all 
1201 If ) corresponds to a stellar density of 6pQ , which 
is inconsistent with main-sequence stellar evolution 
for all stars but late M-dwarfs. This implausibly 
high density derived from a circular orbital fit to 
the light curve implies that the planet has an ec- 
centric orbit and is transiting near periapse (e.g. 
Figure 1 of Paper I) . 

3. The in-transit data feature structures that may be 
caused by star spot crossings (e.g. the bump in the 
purple, solid circle light curve marked C=377.739 
in Figure [1]). The ratio of scatter inside of tran- 
sits to that outside of transits is about 1.2. If 
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Fig. 1. — Detrended light curves, color-coded by transit epoch, 
spaced with arbitrary vertical offsets. The top eight light curves 
are phased based on a constant, linear ephemeris (Table[2] column 
3), revealing the large TTVs. Each light curve is labeled 'C with 
its best-fit mid-transit time (Table [2] column 3). In the second- 
from-the-bottom compilation, each light curve is shifted to have an 
individual best-fitting mid-transit time at t=0. The bottom points 
are the residuals multiplied by 10. Solid lines: best-fitting eccentric 
model (Table [2] column 3). 

the rotation of its spot pattern, we may be able 
to measure the stellar rotation period and com- 
bine it with other stellar parameters to constrain 
the line-of-sight compo nent of the system's spin- 
orbit configuration (e . g. ISanchis-Ojeda et alj|2011t 
INutzman et al.l [20111 IDesert et al ' 
star's surface temperature were 



Mil". 



If the 
gr eater than or 

equal to the KIC estimate of 6498 K (jBatalha et al.l 
I2010D , we might expect the star to l ack a convective 
envelope (|Pinsonneault et al"1l2001l ) and star spots. 
Therefore the star may be significantly cooler than 
this estimate 

The light curve implies that the transiting candidate 
KOI-1474.01 may be an eccentric planet experiencing 
perturbations from an unseen companion and with a 
measurable spin-orbit alignment, an ideal testbed for the- 
ories of planetary migration. However, in order to val- 
idate and characterize the candidate, first we must pin 
down the stellar properties and assess the probability 
that the apparent planetary signal is a false positive. 

3. HOST KOI-1474, A RAPIDLY-ROTATING F STAR 

The stellar properties of KOI-1474 are essential for val- 
idating and characterizing the transiting candidate, but 
the parameters in the KIC are based on broadband pho- 
tometry and may b e systematically in error, as noted by 
I Brown et al.l ( 20111 ). Here we use a combination of spec- 
troscopy (i )3.1[) . photometry ( £|3.2[) . and stellar evolution 
models ( §3.3|) to characterize host star KOI-1474. 

3.1. Stellar temperature, metallicity, and surface gravity 
from spectroscopy 

We obtained two high signal-to-noise, high resolu- 
tion spectra for KOI-1474 using the High Resolution 
Echelle Spectrom eter (HIRES) on the Keck I Telescope 
ijVogt et al.lll994D . The spectra were observed using the 
standard setup of the California Planet Survey, with the 
red cross disperser and the 0'.'86 C 2 decker, but with th e 
iodine cell out of the light path ((Johnson et al.ll2012[ ). 
The first observation was made with an exposure time of 
270 seconds, resulting in a signal-to-noise ratio (SNR) of 
ss 45 at 6000 A; the second exposure was 1200 seconds 
long, resulting in a SNR ss 90. 

As described in Paper I, we use SpecMatch to compare 
the two spectra to the California Planet Survey's vast 
library of spectra for stars w ith parameters from Spec- 
troscopy Made Easy ( SME; iValenti fc Piskunovl 119961: 
iValenti fc Fisch er 2005) . The closest- matching spectrum 
is that of HD3861. In order to match KOI-1474 to this 

9 However, Hirano et al. (2012) recently found photometric vari- 
ability due to star spots for several hot stars, including KOI-1464, 
which has a surface temperature of 6578 ± 70 K, so the signatures 
of star spots we notice are not necessarily inconsistent with KOI- 
1474's KIC temperature. 
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relatively slowly rotating F dwarf, we must rotation- 
ally broaden the spectrum of HD3861. The total line 
broadening for KOI-1474, u rot sini s = 13. 6± 0.5 km/s, 
is a combination of the HIRES instrumental profile, ro- 
tational broadening, and broadening due to turbulence 
(macroturbulence being the dominant term , rather than 
microturbulence: [Valenti fe Fischeil I2005D . We assume 
that KOI-1474 has the same macroturbulent broadening 
and instrumental profile as HD 3861. Then we apply ad- 
ditional rotational br oadening to HD 3861 using MORPH 
([Johnson et al.l 120061 ) to match the spectra of KOI -1474 
using the rotational broadening kernel described bv lGrayl 
(2008]). The v rot sini s for KOI-1474 is 



V T0 t sin i s 



J (v ro t smi s ) 



2 

HD3861 



(v Tot sini s )2 road 



where (v TO t sini s )HD386i = 2.67 km/s is the k nown 
v mt sini s of HD 3861 ([Valenti fe Fischerl [2001 and 
(i> ro t sini s )broad = 13.3 km/s is the additional rotational 
broadening applied to the HD 3861 spectrum to match 
the lines of KOI-1474. See lAlbrecht et al.l ([20ll . §3.1 
for a discussion and demonstration of this technique for 
measuring v TO t sini s . 

Next, from a weighted average of the properties of 
HD 3861 and the other best match spectra, we mea- 
sure an effective temperature T e g = 6240 ± 100 K, sur- 
face gravity logg = 4.16 ± 0.20, and iron abundance 
[Fe/H] = 0.09 ± 0.15. These measured values are con- 
sistent with the KIC estimates of T off = 6498 ± 200 K 
an d log q = 4.08 ± 0.4 (with uncertainties estimated 
bv iBrown etldl 120111 ) but are more accurate and pre- 
cise because they come from high-resolution spectroscopy 
rather than broadband photometry. Based on the re- 
vised, cooler value for its effective temperature, KOI- 
1474 may indeed have a convective envelope and thus the 
structures in the transit light curves (Figure [J) could be 
due to spots. Therefore spot-induced photometric vari- 
ability may allow us to measure the stellar rotation pe- 
riod Prot f ^3.2[) . which we can combine with other stellar 
parameters to infer the transiting candidate's projected 
spin-orbit alignment ( EI5.2j) . 

3.2. Stellar rotation period from photometry 

KOI-1474 appears to exhibit rotational photometric 
variability due to star spots, which cause the star to ap- 
pear brighter (dimmer) as the less (more) spotted hemi- 
sphere rotates into view. We see what may be an effect of 
these spots in the purple, solid circle light curve marked 
C=377.739 in Figure [TJ a bump during transit consistent 
with a planet crossing a star spot. A periodogram (not 
shown) of the entire photometric dataset (Q0-Q6) ex- 
hibits a prominent cluster of peaks near 5 days. However, 
a periodogram is not the best tool to measure stellar ro- 
tation periods because: a) the photometric variability is 
non-sinusoidal, and b) the spot pattern is not expected to 
remain coherent over the entire 508-day dataset and thus 
the phase and amplitude of the best-fit sinusoid change 
over the data's timespan. 

To obtain an optimal measurement of the stellar rota- 
tion period, we comp ute the discrete-correlat ion-function 
(DCF, Equation 2 of lEdelson fe Krolikj[l988l) . which was 
rec ently used to meas ure the rotation p eriod of Corot- 
7 (|Queloz et al.l I2009D and Kepler-30 ([Fabrvckv et al.l 
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Fig. 2. — Discrete-corrclation-function (DCF, Edelson and Krolik 
1988) for the long-cadence Kepler Q0-Q6 KOI-1474 photometric, 
dataset as a function of time lag. The peak at 4.6 ± 0.4 days 
corresponds to the stellar rotation period. 
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Fig. 3. — Posteriors (solid) of stellar radius (panel 1), mass (panel 
2), and density (panel 3) in solar units. The posteriors obtained 
from the prior alone (dashed gray) and from the data alone (dotted) 
are plotted in each panel, demonstrating that our data provide 
stronger constraints on the stellar parameters than do our priors. 



2012). First we detrended the data with the PyKE rou- 



tincdj using co-trending vectors. iWelshl (j!999f ) found 
that it is crucial to remove long-term trends from the 
time series before applying the DCF or biases may re- 
sult. Then we computed the DCF using the Institut 
fur Astronomie und Astrophysik Tubingen DCF rou- 
ti neF 1 ] an IDL implement ation of the DCF described 
in lEdelson fc Krolikl ([1988) . The possible range for the 
DCF is -1 to 1; the amplitude is normalized such that 
DCF = 1 indicates perfect correlation. We plot the DCF 
(computed with a lag range of 0.1 days to 20 days and 
with 200 frequencies) as a function of time lag in Figure 
[H The DCF is highest in the region lag < 0.2 days (i.e. 
lags that are small but greater than 0, for which the DCF 
=1 by definition), indicating that most of the photomet- 
ric variability occurs on short timescales, most likely due 
to a combination of high-frequency stellar variability and 
instrumental noise. However, we also see lower amplitude 
but pronounced peaks at longer periods. 



Guest 



10 Available at K epler 
http: / /keplergo. arc. nasa.gov 

w Available at 

tuebingen.de / software / idl / aitlib / index. shtml 
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The DCF exhibits the variations we expect due to star spots. Imagine observing the star at time t; the hemisphere 
in view has either more or fewer spots than the unseen hemisphere. At time t + P rot /2, the other hemisphere has 
fully rotated into view, so the flux at t and t + P rot /2 are negatively correlated. Therefore, we interpret the negative 
DCF near 2 days as corresponding to half the stellar rotation period. At time t + P ro t, we see the same hemisphere 
as at time t; therefore, we interpret the strong positive correlation at lag 4.6 ±0.4 days as the stellar rotation period, 
for which uncertainty range corresponds to the width at half-maximum. The amplitude DCF = 0.1 indicates a 10% 
correlation between points separated in time by P rot . The other hemisphere rotates fully into view again at i+3P ro t/2, 
corresponding to the negative DCF at lag 7 days; at lag 2P rot = 9 days, the DCF is positive again. This pattern 
continues, and the amplitude would remain constant if the spot pattern were constant. However, the spot pattern is 
changing over time, so the amplitude of the correlation "envelope" decreases with ti me lagE] The measured ro tation 
period of 4.6 days is consistent with the distribution measured for F, G, K stars bv lReiners fc Schmittl (|2003[) ; they 
find that the distribution of projected rotation periods (i.e. the rotation periods measured from i> rot sin i s assuming 
i s = 90°) peaks at 5 days. 

3.3. Stellar density from evolution models 

The candidate's orbital eccentricity, the ultimate quantity of interest, depends weakly on the host star's density 
(see Paper I and references therein). Thus it is important to have an accurate, if not precise, estimate of the host 
star's density and, importantly, a conservati ve estimate of the un certainty. For this task we use the finely-sampled 
YREC stellar evolution models computed bv iTakeda et al.l (|2007t ). sampled evenly in intervals of 0.02 dex, 0.02 Mq, 
and 0.02 Gyr for metallicity [Fe/H], stellar mass M+ and age r* respectively. The model parameters are stellar age 
t*, mass M^, and fractional metallicity Z, and we wish to match the effective temperature T e g, surface gravity logg, 
and [Fe/H] measured spectroscopically in £13.31 along with their 68.2% confidence ranges denoted by their "one- 
sigma errors" {or oft , o-\ ogg , cr\p e /m}, respectively. In what follows, the subscript "spec" refers to the spectroscopically 
measured quantity, while quantities with no subscript are the model parameters. 

Applying Bayes' theorem, the model posterior probability distribution is 

prob(M*, r*, Z|T c ff !Spoc , [Fe/H] spcc , logg spoc , I) oc 
prob(T eff!Spcc , [Fe/H] spcc , log 5spec |M*, r*, Z, I)prob(M*, r*, Z\I) (1) 

where / represents additional information available to us based on prior knowledge of the Galactic stellar population. 

The first term on the right hand side (RHS) is the likelihood, which we compute by comparing the effective 
temperature, surface gravity, and metallicity generated by the model to the values we measured from spectroscopy: 



prob(T c ^ spcc [Fe/H] spcc ,log 3spcc |M„r„Z,7) cx exp [-% j exp ( ) exp [ j (2) 



where 



-v 2 

X[Fe/H] 



2 [T(M^, Tj,, Z) - Teff.spcc] 

XT eff - ~2 

u Toff, spec 

[Fe/H](M„r„Z)-[Fe/H] spcc - 2 



2 

CT [Fo/H],spec 

[log g(M+,T„Z) - log 

Xlogg — ~2 [3) 

log g.spcc 

The second term on the RHS of Equation fT]). prob(iiP, r*, Z\I), is the prior information known about the 
model parameters. Here we make use of some additional information I — the galactic latitude and longitude of 
the Kepler field and the measured apparent Kepler magnitude of KOI- 1474 — to infer the relative probability of 
observing different types of stars. A number of factors go into this probability, including the present-day stellar mass 
function, the volume distribution and ages of stars along our line of sight to the Kepler field, and t he Malmquist 
bias. Fortunately, the TRILEGAL code (TRIdimensional modeL of thE GALaxy; iGirardi et al.l[20 05) synthesizes a 
large body of observational, empirical, and theoretical studies to produce a model population of stars in the Kepler 
field that are consistent with KOI-1474's apparent Kepler magnitude Kp = 13.005 ± 0.030 (Batal ha et al.l [2010) 
and Galactic coordinates. From this model population, we use a Gaussian kernel density estimator to compute a 
three-dimensional density f unction for the prior prob(M*, t*, Z\I). 

Each combination of Takcda et al.l (|2007t ) model parameters — (M+, Ti,,Z) — has a corresponding P* and L*, 

and we calculate the corresponding stellar density = jg£ (-^) 3 Po- We compute the star's absolute Kepler 
bandpass magnitude -?Cp,absoiute through the follow steps: we transform L± into a V magnitude using a bolometric 
correction, transform V into the absolute Sloan magnitude g, and compute the distance modulus using the difference 

Unfortunately, the d ecrease in the correla t ion am plit ude with lag imp l ies th at we are unlikely to be able to measure the stellar obliquity 
a using the method of Sanchis-Oicda ct al. (2011) and Nutzman ct al. (2011). The spot cycle is likely shorter than the interval between 
subsequent transits. 
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between the absolute g magnitude and the a pparent g 
magnitude from the KIC (jBatalha et al.ll2010D . Then we 
apply the distance modulus to the apparent Kp to ob- 
tain -ftTp, absolute- Thus we can transform the model pos- 
terior prob(M*, r*, Z\T c s ySpcc , [Fe/H] spcc , logg spoc , I) into 

posteriors for the stellar stellar properties M*, r*, i?*, p*, 
L*, and ^p, a bsolute (Tabled! column 3). In Figure |3] we 
plot the resulting posteriors for M*, R±, and p*. We 
also plot the same distributions obtained from the data 
alone and from the priort(3 alone; evidently most of the 
constraint comes from the data (i.e. the spectroscopic 
quantities) . 

The derived density for KOI-1474, 0.44+^/Jq, has an 
unc ertainty range encom passing the KIC value of 0.26 
Po (jBatalha et al.ll2010D . The star is significantly less 
dense than the valu e of &pp, derived from a/R+ in the 
table of candidates ()Borucki et al.l l201~li : Batalh a~et al.l 
120121) . Therefore, planet candidate KOI-1474.01 is likely 
to have a large eccentricity, which we will measure in 
SJS] Fortunately, as shown as Paper I, even the loose 
constraint on the stellar density derived here will result 
in a precise measurement of the candidate's large orbital 
eccentricity. 




0.5 1.0 1.5 2.0 



angular separation (arcseconds) 

Fig. 4. — Sensitivity to off-axis sources in the immediate vicin- 
ity of KOf-1474 using adaptive optics imaging observations with 
NfRC2 at Keck in the K'-band (A c = 2.12/xm). 



4. FALSE POSfTfVE PROBABILITY 

Although a transiting planet may cause the photomet- 
ric signal observed in light curves (Figured]), any of sev- 
eral scenarios involving stellar eclipsing binaries might 
cause a similar signal. This is the well-known problem 
of astrophysical false positives for transit surveys (e.g. 
IBrown|[2003l : iTorres et al.H201lD . Traditionally transiting 
planets have been confirmed through detection of their 
radial velocity (RV) signals. However, the Kepler mission 

12 The Mi, prior probability appear s truncated be l ow M t = 0.78 
in Figure l3l because we only compute Takcda et al. (2007) models 
above this value. However, the likelihood completely rules out stars 
with < IMq . 




log(a) r[days] 



Fig. 5. — Three-dimensional probability distribution for the 
trapezoidal shape parameters (depth <5, duration T, and "slope" 
T/t) for the nominal planet scenario. The distributions are gen- 
erated by sim ulating a statistically representative population (see 
IMortonll2012j . §3.1) for the scenario and fitting the shape param- 
eters to each simulated instance. Each population begins with 
100,000 simulated instances, and only instances that pass all avail- 
able observational constraints are included in these distributions. 
In this case, no additional observational constraints are available so 
the 100% of the distribution remains. The transit's shape parame- 
ters S, T, and T/t are marked on each plot with an "X" denoting 
the the median of an MCMC fit. 




\og[S) r[days] 



Fig. 6. — Same as Figure |51 for the HEB scenario. In this case, 
the upper-limit of 200 ppm we place on the secondary eclipse depth 
eliminates 30.4% of the distribution and limits from the Ks-band 
adaptive optics image eliminate 29.1% of the distribution, leaving 
40.5% remaining. 

has necessitated a different paradigm: one of probabilis- 
tic validation. If the false positive probability (FPP) of 
a given transit signal can be shown to be sufficiently low 
(e.g. <C 1%), then the planet can be considered validated, 
even if not dynamically confirmed. Here we attempt to 
validate KOI-1474.01 but find a 3.1% probability that 
the signal is due to an astrophysical false positive. 

At first glance, the short duration of KOI-1474.01's 
transit (Section [5]) causes particular concern: the sig- 
nal could be a transit or eclipse of an object orbiting a 
smaller, blended star, which would make the duration 
more in line with that expected for a circular orbit. In 
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TABLE 1 

Stellar Parameters for KOI 1474 



Parameter 


Value a 






Measured 


Derived from model 


Right ascension, RA (hour,J2000) 


19.694530 




Declination, Dec (degree, J2000) 


51.184800 










Stellar effective temperature, T fr [K] 
i en L J 


6240±100 


6230±100 


Iron abundance, [Fc/H] 


0.09 ±0.15 


00 + 16 
uuu -0.12 


Surface gravity, log(^[cms ^1 


4 16+0 20 


4 23+ - 13 


Limb darkening coefficient, 




0.320 ± 0.015 


Limb darkening coefficient, /^2 b 




0.304 ± 0.007 


Main sequence age, t* [Gyr] c 




2 8+H 
1.2 


Stellar mass, M* [M0] c 




-, 99 +0.12 

^--^-o.os 


Stellar radius, [-Rq] 




X.tU_ Q 21 


Stellar density, P* \pa\\ 




0.44+O^f, 


Stellar luminosity, L* [-Lq] 




9 7+1-6 
z - '-0.8 


Apparent Kepler-band magnitude, Kp 


13.005 ± 0.030 


R +0.4 
d -°-0.5 


Absolute Kepler-band magnitude, -ft"p, a bsolute 




Distance (kpc) 




0.78l°; 2 3 


Rotation period, P rot [days] 


4.6 ± 0.4 


14.7i?:« 


Rotation speed, w r ot[km s _1 ] 




Sine of stellar spin axis inclination angle, sin i s 




0.93+^4 


Stellar spin axis inclination angle, i s [degree] 






Deviation of stellar spin axis from edge-on, |90 — is [degree] 




21±il 



a The uncertainties represent the 68.3% confidence interval of the posterior distribution. 
b Sing 2010 



A prior was imposed on this parameter. 

KOI 1474.01 

Priors Likelihoods 



eb: 1.4e-03 
heb: 1.2e-03 
bgeb: 3.6e-05 
bgpl: 1.5e-04 
pi: 2.7e-04 




eb: 3.3e-04 
heb: 6.0e-03 
bgeb: 3.2e-ll 
bgpl: 5.4e-14 
pi: 8.7e 01 





eb: 0.002 




heb: 0.029 




bgeb: 0.000 




bgpl: 0.000 




pi: 0.969 


fpiy 


=0.064 



Final Probability 

Fig. 7. — Prior (top left), likelihood (top right), and final (bot- 
tom) probabilities for four false positive scenarios — an undiluted 
eclipsing binary ("eb"), hierarchical eclipsing binary ("heb"), back- 
ground eclipsing binary ("bgeb"), and background planet ("bgpl"). 
The priors and likelihoods are computed following Morton (2012). 
Each final probability is the product of the scenario's prior and like- 
lihood, normalized so that the total probabilities sum to 1. The 
quantity f p i^y indicates the specific occurrence rate for planets of 
this size that we would need to assume in order for the FPP to be 
less than 0.5%. Since this rate, / p i,v = 6.4%, is higher than our 
assumed / p i = 0.01, we do not consider the candidate validated. 



order to calculate the F PP for KOI-14 74.01, we follow 
the procedure outlined in iMortonl (|2012t ). which incorpo- 
rates simulations of realistic populations of false positive 
scenarios, the KIC colors, the measured spectroscopic 



stellar properties, and a descriptive, trapezoidal fit to 
the photometric data. 

To place constraints on blending by searching for 
nearby sources, we obtained adaptive optics images of 
KOI-1474 on March 29, 2012 using NIRC2 (PI: Keith 
Matthews) at the 10m Keck II telescopes. KOI-1474 
is sufficiently bright to serve as its own natural guide 
star (Kp = 13.005) and therefore does not require the 
use of a laser to correct for wavefront errors introduced 
by the Earth's atmosphere. Our observations consist of 
18 dithered images (10 coadds per frame, 2 seconds per 
coadd) taken in the K' filter (A c = 2.12/im). We used 
NIRC2's narrow camera mode, which has a platescale 
of 10 mas / pix, to provide fine spatial sampling of the 
stellar point-spread function. 

Raw frames were processed by cleaning hot pixels, flat- 
fielding, subtracting background noise from the sky and 
instrument optics, and coadding the results. No off- 
axis sources were noticed in individual frames or the fi- 
nal processed image. Figure [4] shows the contrast levels 
achieved from our observations. Our diffraction-limited 
images rule out the presence of contaminants down to 
AK ' = 5 mag and AK' — 8 mag fainter than the primary 
star for separations beyond 0.2" and 0.7" respectively. 

We plot the probability distributions for the nom- 
inal planet scenario in Figure [5l as well as for the 
most likely alternative to a transiting planet: an 
hierarchical eclipsing binary (HEB) (Figure in 
which KOI-1474 has a wide binary companion of 
comparable brightness (within a few magnitudes) 
that is being eclipsed by a small tertiary companion. 
The probability of the HEB scenario is 2.3%. In Figure[71 
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TABLE 2 

Planet Parameters for KOI 1474.01 



Parameter 


Value a 






Circular fit 


Eccentric fit 



Average mid transit epoch, 
Mid transit epoch of transit 1, 



Average orbital period, P [days] 



days 




BJD-2455000 


days 




BJD-2455000 


days 




BJD-2455000 


days 




BJD-2455000 


days 




BJD-2455000 


days 




BJD-2455000 


days 




BJD-2455000 


days 




BJD-2455000 


days 




BJD-2455000 



T c 
Tl 
T-2 
T 3 
Ti 
T 5 
T 6 
T 7 
T a 

Planet- to-star radius ratio, Rp/R t 
Stellar density, p* 
Orbital inclination, i [degree] 
Limb darkening coefficient, fi± c 
Limb darkening coefficient, fi2 c 
Impact parameter, b 
Planetary radius, R p [i?®] 
Normalized red noise, o> 
Normalized white noise, tj^j 
Eccentricity, e 

Orbital period after tidal circularization, Pfl na i 
Line-of-sight spin-orbit angle, \i — i s \ [degree] 



69.7339±0.0016 

238.273±0.0I1 

-40.6701±0.0008 

29.0600±0.0006 

98.7647±0.0006 

168.5752±0.0006 

238.3146±0.0005 

308.0092±0.0008 

377.7250±0.0006 

447.4555±0.0006 

n rims +0 0007 
0.U618 _ 0003 

+0.4 
1.6 
+0.05 
-0.08 

n -J14+0 018 
u.ji4_ 012 

302 + 006 
+0.21 
0.12 



9.2 
'.93 



0.18] 



00005+"' "" 07 

U.UUUUO_ odoi);) 

o nonm +0 - l)()l)()1() 

U.UUU1J1_ .000004 



69.7340±0.0015 
238.273±0.010 
-40.6702±0.0009 
29.0600±0.0007 
98.7647±0.0007 
168. 5752±0. 0007 
238. 3146±0. 0007 
308. 0092±0. 0009 
377.7250±0.0007 
447.4555±0.0007 

+0.0006 
-0.0004 
K c+0.30 
-0.10 
+0.4 
-1.3 

U.Oll_ 012 
f0.005 
-0.009 
4+0.25 
-0.09 
-+2.4 

00007+ - 00005 

U-UUUU/ _ . ()()()().-, 
,,+0.000007 
-0.000007 

n si +010 

0.»l_o.07 
/I+6 
-10 

-17 
-14 



0.0617 
0.36 c 
89.2 
lj 

0.304j 
0.14^ 
9.5] 



0.000134] 



11] 
21. 



a The uncertainties represent the 68.3% confidence interval of the posterior distribution. 

b P and T c are determined from a linear fit to the transit times. The uncertainty in T c is the 

median absolute deviation of the transit times from this ephemeris; the uncertainty P is this 

quantity divided by the number of orbits between the first and last observed transits. 

c A prior was imposed on this parameter. 

summarize the prior, likelihood, and total probability of the nominal transiting planet scenario compared to that of 
several false positive scenarios. The FPP is: 



FPP 



L 



FP 



L 



FP 



(0.002 + 0.029 + 0.000 + 0.000) 

Min 969 



jf£r L T p (0.002 + 0.029 + 0.000 + 0.000) 



= 0.031 



(4) 



where Lpp is the sum of the probabilities of the false-positive scenarios, Lpp is the probability of the nominal planet 
scenario, and fp is the assumed specific occurrence rat^B for planets between 5.7 and 11.3 i?®. Although this FPP 
is low, we do not consider it sufficiently low to validate the planet. In the analysis following in the remainder of the 
paper, we assume that KOI-1474.01 is a planet and refer to it as "planet," but in fact it remains a candidate planet. 
We are conducting a radial-velocity follow-up campaign of this target to confirm this candidate by measuring its 
mass. 



5. THE HIGHLY ECCENTRIC ORBIT OF KOI-1474.01 

In <JH we revised the stellar properties of KOI-1474 and found that the star's density indicates that the 
we(validated) planet's orbit is highly eccentric. To qu antify the eccentric ity, we now model the light curves (Figure 
[1| with the Transit Analysis Package software (TAP, G azak et al.ll2012T ) to obtain the posterior distribution for the 
eccentricity and other transit parameters ($5T}, using the technique described in Paper I. In £15.21 we place constraints 
on the spin-orbit alignment based on stellar properties measured in £13.31 In H5.31 we assess the observed TTVs and 
explore the nature of the third-body perturber. 



5.1. Fitting orbital parameters to the light curve 

Here we measure KOI-14 74. Ol's orbi t al par ameters, including eccentricity, from the transit light curves (Figure 
J). We use TA P to fi t a IMandel fc Ago! (I2002T) light curve mo del, employ ing the wavelet likelihood function of 
(barter fc Winn! ((2009). We replace the parameter a/R* with (|Winnll2010l Equation 30) in the limit that (ill* + 



M p )/(|7ri?*) 3 — > p+, but transform p+ into a/ to compute the light curve model. Using the spectroscopic stellar 
parameters measured in H3.ll (Table HI column 2), we calculate th e limb dark ening coefficients pi and /i 2 and their 
uncertainties with the table and interpolation routine provided by Sing (201(|). In all the orbital fits discussed 

The assumed 1% occurrence rate is motivated by the debiased 1% occurrence rate for hot Jupit e rs in t he RV sample (Wright et alJ l2012T ). 
a In order to produce a FPP of than 0.5%, f p would have be greater than 6.4%. See [Morton (2012) for a discussion of specific planet 
occurrence rates. 
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here, we impose normal priors on pi and \ii (Table [T]), which are well-measured for the Kepler bandpass. We also 
verified that uniform priors on the limb darkening yield consistent results (with slightly larger uncertainties) for all 
the orbital fits we perform. The other light curve parameters we fit for are the mid-transit time of each light curve 
T, the planet-to-stellar radius ratio i? p /-R*, the fractional white noise a w , the red noise ay, the inclination i, and the 
argument of periapse u>, with uniform priors on each of these quantities. 

Finally, to speed up the fit convergence, we explore parameter space using the parameter g instead of the 
planet's orbital eccentricity e. The parameter g corresponds approximately to the ratio of the observed transit speed 
to the speed expected of a planet with the same period but e = 0: 

ii- / \ V 3 

1 + e sin lu x 



--circ 



We impose a prior on g to maintain a uniform eccentricity prior (see §3.3.1 of Paper I for further details): 



sin 2 lu (sin 2 lu — l) + g 2 (l + sin 2 lu) ± 2g sin lu y sin 2 lu — 1 + g 2 
P rob (f ) = ^T7~ ——2 — (6) 



\J sin 2 lu - 1 + g 2 (g 2 



sm 2 lu I 



for which the + corresponds to g > 1 and the — to g < 1. We transform g into e to compute the light curve model. 

First we fit a circular orbit (Table [3J column 2), fixing e = and leaving free p*, to which we refer as 
Pcirc- We find that: 1) although we only have long-cadence data for KOI-1474 (Figure [T]), p circ and the impact 
parameter b are separately well-constrained (see also §4.2 of Paper I for a discussion of long-cadence data), and 
2)the p* posterior computed from stellar properties in 33.31 (p+ = O.44t^2o/ 9 0) ^ aus ^ ar outside the transit light curve 
posterior distribution for p circ (p* = 9.2t?;gPO' Figured! top left panel), where the uncertainties indicate the 68.3% 
confidence interval. Thus a circular fit is inconsistent with our prior knowledge of the stellar parameters. 

Because the eccentricity depends only weakly on the assumed stellar density, the eccentricity measurement we 
are about to perform is relatively robust to errors in the assumed stellar density. When p* > p c irc, the transiting 
planet has a minimum eccentricity obtained by setting lu = ir/2 in Equation ^ (i.e. the planet transits at periapse). 
Imagine that p* were biased or in error. The fractional change in e m i n would be: 



Ae ■ 4 M—) 



/ ^ x2/3 / ^ \-2/3 

\ Pcirc J \ Pcirc J 



(_£*_) 

\ Pcirc > 



(7) 



The ratio = = 21, corresponding to e min = 0.77 and ^ emi " = 0.18 /pS i c V . So if the stellar density were 

^circ 

biased upward by 10%, the minimum eccentricity would be biased upward by only 1.8%. See §3.1 and §4.1 Paper 
I for a detailed exploration of how the stellar density's assumed probability distribution affects the eccentricity 
measurement. 

Next we fit the light curve allowing the planet to have an eccentric orbit (Tabled column 3) and using the 
stellar density posterior from £13.31 as the stellar density prior for the light curve fit. As argued in Paper I (Section 3), 
an MCMC exploration — as implemented in TAP — naturally accounts for the transit probability and marginalizes 
over the uncertainties in other parameters. Even though e and lu are degenerate for a given g (Equation [5]) , there is 
a lower limit on e, and the posterior falls off gradually, as e — > 1 and the range of possible lu satisfying Equation [5] 
narrows. The posterior distributions for e and lu are plotted in FigureO We measure e = 0.81^q g?- F° r comparison, 
if we had set the stellar density prior to be uniform between O.lp© — O.2/?0 (0.6p© — 1.2p0), we would measure 
e = 0.90t ;° 3 3 (e = 0.73±° : Jg). 

By conservation of angular momentum, this planet would attain a final period Pfi na i(l — e 2 ) 3 / 2 = 14ti days if 
it were to undergo full tidal circularization. In §33 we will discuss whether the planet is best classified as a proto-hot 
Jupiter — likely to circularize over the star's lifetime and achieve a short-period orbit — or a failed-hot Jupiter, just 
outside the reach of fast tidal circularization. 



5.2. Constraints on spin-orbit alignment 

Whatever process perturbed KOI-1474.01 onto an eccentric orbit may have also tilted the planet's orbit from 
the plane in which it formed. With a temperature of 6240 ± 100 K (Section GO}, KOI-1474 sit s right on the 625 
K boundary between hot stars with high obliquities and cool stars with well-aligned planets ()Winn "et^l20Tfi) . 
However, if 1) cool st ars have low obliqu ities because their hot Jupiters have realigned the star's outer convective 
layer, as proposed by IWinn et all (|2010D . and 2) KOI-1474.01 is a failed-hot Jupiter, with a tidal dissipation rate 
too low to experience significant circularization over KOI-1474's lifetime, then KOI-1474.01 may have also not yet 
realigned KOI-1474's outer layer. 
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Ultimately we will wish to determine tp, t he total misalignment between the orbit normal and the host star spin 
axis, from three measured projected angles (jFabrvckv fc Winnl l2009; Schlaufman 2010): i, the inclination between the 
planet's orbit and the observer's line of sight; the sky-projected spin-orbit angle A; and i s , the inclination between the 
stellar spin axis and the line of sight. We measured i from the transit light curve in §5.11 (Tabic [21) . The sky-projected 
spin-orbit angle A could one day be measured via the Rossiter-McLaughlin (RM) effect (|McLaughlinlll924t IRossiterl 
Il924t iQueloz et al.l l2000h . the change in the observed radial velocity as a transiting planet blocks por tions of the 
star rotating toward or away from the observer. The effect has a maximum amplitude of about 50 m/s ()Winnl [20101 
Equation 40) However, because KOI-1474.01's transits can occur early or late by over an hour, RM measurements 
of KOI-1474.01 will remain challenging until the TTV pattern "turns over" in future Kepler observations, allowing 
us to predict future transits to much higher precision (Section I5.3j) . We can measure the third projected angle, i s , 
from (u rot sint s ) spec (Section |3.1[) and the posteriors of P rot (Section |3.2[) and P* (Section 13.31) . an approach that was 



recently applied bv lHirano et al.l <|2012f ) to fifteen KOI systems. KOI-1474's rotational velocity is w ro t = 2 p R * and we 
have measured the projected rotational velocity (u ro t sini s ) spcc . Therefore we can find the angle of that projection, 
i s . According to Bayes theorem: 



prob ((vrot sini s ) spoc |P rot , P*, i s ) prob (P rot , P*, i a ) . 



(8) 



prob(P ro t,P*,is|( ^rot sini s ) S p ec ) 

The prior, prob (P rot , P*, « s ), is the product 

prob (P rot , P*, i s ) = prob(P rot )prob(P*)prob(i s ) 

where prob(P rot ) is a normal distribution with mean 4.6 m/s and standard deviation 0.4 m/s ( £13.21) and prob(P*) is 
the posterior from H3.3I Assuming stellar spin axes are randomly oriented throughout the Galaxy, the distribution 



of cosi s is uniform and thus prob(i s ) = £ sm! s . 

Next, we integrate Equation [8] over P ro t and P* to obtain the stellar inclination i s conditioned on our measured 
projected rotational velocity v rot sini s . 



prob (i s 



sim 



;)spec) 



prob ((v Iot sinj s ) spec |P rot , P*, i s ) prob (P rot , P*, i s ) dP rot dR+ 



(9) 



As a practical implementation of Equation ^ we randomly draw P rot and P* from the distributions calculated in 
^3.21 and M3 .31 respectively and i s from a uniform distribution of cosi between and 1. Drawing from these respective 
distributions is equivalent to creating a grid in these parameters and subsequently downsampling according to the 
prior probabilitities. Then we compute the likelihood 



prob((u r otsin« s ) S pec|Prot,P*,« s ) = cxp 



27rP* 



■ am. 



(v m t sini s ) s 



2a 



(t> rot sini„) s 



(10) 



where (v 

rot sin i s )spcc — 13.6 m/s and &(y TO x sinz s ) spcc — 0.5 m/s (Section 13. ip . Then we select a uniform random 
number between and 1; if the uniform random number is less than prob ((i> ro t sinz s ) spec |Prot, P*, is) (Equation [TO)), 
we include the model (P ro ti P* ; *s) m the posterior. We repeat drawing (P ro t, P*, i s ) models until we have thousands 
of models that comprise the posterior. 

We measure a projected angle for the stellar spin axis i s = 69l^ 7 degrees. Combining the posterior of i s with 

the posterior of the planet's inclination i (Section 15. II) .we obtain \i — i s \ = 21~^\\, for which the total uncertainty is 
dominated by the uncertainty in the stellar radius. We list these angles in Table [U and plot the posterior for the 
line-of-sight spin-orbit angle \i — i s \ in Figure [8] (top right panel). Our posterior distribution is consistent (within 
2 a) with close alignment, yet allows misaligned configurations as well. We also caution that differe ntial rotation 
may cause systematic errors in the measured alignment, depending on the latitude of the spots (see IHirano et all 
2012, Section 5.3 for a detailed discussion). Furthermore, the line-of-sight spin-orbit angle |z — i s \ offers no constraint 
on whether the planet's orbit is prograde or retrograde. However, two types of future follow-up observations will 
allow us to better constrain the planet's orbit in three dimensions. First, additional constraints on the planet's orbit 
through radial- velocity measurements will in turn constrain the stellar radius, providing a more precise measurement 
of \i — i s \. To this end we are currently conducting a Doppler follow-up program at Keck with HIRES. Second, 
from the measurement of the sky-projected spin-orbit angle A via the RM effect, the total spin-orbit angle "0 can be 
computed by combining A with a refined line-of-sight measurement \i — i s \. 



5.3. Transit Timing Variations 

The light curves in Figure Q] reveal large variations in the mid-transit times of KOI-1474.01, which may 
be caused by perturbations from another planet or sub-stellar companion. If KOI-1474.01 underwent HEM, this 
perturber may have been responsible. Table [2] displays the mid-transit times from the orbital fits performed in 
Section 15.11 There the best-fitting linear ephemeris is also given, from which the times deviate significantly. In 
Figure O we plot an observed minus calculated (O-C) diagram of the observed transit time minus the transit time 
calculated from a constant orbital period. The scale and sharpness of the features in Figure [S] suggest a nearby giant 
planet or brown dwarf perturber. We assume this perturber is on an exterior orbit, as KOI-1474.01's eccentric orbit 
leaves little dynamical room interior to itself. We are undertaking a radial-velocity follow-up campaign (Johnson et 
al. 2013, in prep) that may allow us to rule out an interior, Jupiter-mass companion. 
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Fig. 8. — Top left: p c i rc obtained from circu lar f it to the transit 
light curve (solid) and posterior for p* from £13.31 (dashed); since 
the host star is not highly dense (i.e. the two posteriors do not 
overlap), the planet's orbit must be highly eccentric. Top right: 
Posterior for projected spin-orbit alignment from an eccentric fit to 
transit light curve, imposing a prior on p*. Bottom left: Posterior 
distribution u from an eccentric fit to transit light curve, imposing 
a prior on p*. Bottom right: Joint posterior for ui vs. e. The black 
(gray, light gray) contours represent the {68.3, 95, 99}% probability 
density levels (i.e. 68% of the posterior is contained within the 
black contour) Over-plotted as a black-and-white dotted line is 
a histogram of the eccentricity posterior probability distribution 
marginalized over u). 

The "jump" in the O-C diagram likely corre- 
sponds to the periapse passa g e of a n eccentric com- 
panion (IBorkovits et all l2003l l200i lAeol et all [20051: 
iBorkovits et al.1 120111 ). Throughout its orbit, this per- 
turbing companion creates a tidal force on the orbit of 
the transiting planet. If the companion's orbit is exte- 
rior to and within the plane of the transiting planet's, 
the tidal force increases the inner planet's orbital period 
or, equivalently, decrease s the effective ma ss of the cen- 
tral star (see Section 4 of lAgol et al.ll2005l for a detailed 
derivation). The tidal force varies with the distance be- 
tween the perturber and star and is strongest when the 
perturber is at periapse. Therefore, as the perturber ap- 
proaches periapse, the transiting planet's orbital period 
lengthens, causing later and later transit arrival times, 
corresponding to the discontinuity seen in Figure [5] The 
period of the TTV cycle corresponds to the perturb- 
ing planet's orbital period. The amplitude is set by the 
change in the tidal force (a combination of the perturbing 
planet's mass and periapse distance, which is a function 
of the eccentricity and orbital period). The sharpness 
of the O-C depends on the perturber's eccentricity - 
whether the perturbation is the flyby of a companion 
on a highly eccentric orbit or the gradual approach of a 
moderately eccentric companion. The transiting planet's 
orbital eccentricity also subtly affects t he shape of the 
O-C d iagram, as explored in detail by IBorkovits et al.l 
(|201lD . Our Figure [9] has a similar appearance to the 



TTVs produced by IBorkovits et al.l (|201lD 's analytical 
and numerical models of eccentric, hierarchal systems. 

Currently we do not have a long enough TTV base- 
line to uniquely mod el the perturbing companion, as 
INesvornv et al.l ()2012D achieved for the system KOF857. 
Since "jumps" in the O-C diagram correspond to the 
perturber's periapse passage and we have only seen one 
such jump, apparently the current TTVs cover less than 
one orbit of the outer companion. Therefore we cannot 
well constrain the outer body's orbital period. The TTV 
amplitude - set by the tidal force on the transiter's or- 
bit - is well constrained but depends on the perturber's 
mass, orbital period, and eccentricity; therefore we ex- 
pect to find degeneracy among these quantities. Fur- 
thermore, the tidal force on the transiter's orbit depends 
on the mutual inclinations of the bodies. The tide due 
to a polar position for the perturber would decrease the 
transiter's orbital period; averaging over the bodies' po- 
sitions, a very inclined perturber could be more massive 
and yet produce a comparable amplitude perturbation. 

Here we explore a subset of all possible parameters for 
the perturbing planet. With only eight transit times (Ta- 
ble [2]) , we have a great amount of freedom in the fits, but 
it is still of interest whether or not a physical model of 
a perturber can fit these data{3 Thus we proceed with 
direct 3-body fits to the data. We do not expect the 
TTVs to be sensit ive to the mass of the t ransiting planet 
or th e host star (pBorkovits et al.1 120111 : INesvornv et al.1 
120121 ) so we fix M m = lM Jup and M* = 1.22M Q . We 
fix the eccentricity and argument of periapse of KOI- 
1474.01 to various values consistent with the light curve, 
then fit for the period P2, the conjunction epoch Tq^, 
e2 cos UJ2 , e2 sin uj2 , and the mass M2 of the perturb- 
ing body (denoted "2"). The fits are performed via a 
Levenberg-Marquardt algorithm drivi ng a numerical i n- 
tegration that solves for transit times (|Fabrvckvll2010( ). 

Initially we consider coplanar, edge-on orbits. This 
configuration is consistent with the transiting planet, and 
although no transit of the perturbing body has been ob- 
served, it may transit in future data or may be within 
a few degrees of edge-on, which would make little dif- 
ference to the TTVs. We first allow all 5 parameters of 
the outer planet to float freely, finding the best fits at 
each value. We performed two fits (Tabled! rows 1-4), 
one with KOI-1474.01 transiting at periapse and another 
with it transiting at semilatus rectum. Both fits are ac- 
ceptable, so we find that we cannot currently use TTVs 
to distinguish these possibilities. In Figure[9j we plot the 
O-C variations generated by these two models. In both 
cases, the perturber is a giant planet on a moderately ec- 
centric orbit with a roughly Martian orbital period. We 
repeated both these fits with a fixed mass of 100Mj up for 
KOI-1474.01 and found, as expected, that the solutions 
were similar, with only a slightly larger (~ 20%) best-fit 
mass for the perturber. 

Next we perform a fit for which the transiting planet 
and the perturbing body have a 124° mutual inclination, 
a possible outcom e of the secular chaos HEM mechanism 
(jNaoz et al.l 120111 ) . As discussed above, non-coplanar 

13 For example, Ncsvornv ct al] H2012T 1 demonstrated that, as 
expected, they could not find a physically plausible model when 
they scrambled their TTVs. Failure to find an orbital model that 
reproduces the observed TTVs would cast suspicion on our inter- 
pretation that they are the signature of an unseen companion. 
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Fi g. 9. — Left: Observed mid transit times (purple dots) of the eight transits of 1474.01 with subtracted best fit linear cphemcris from 
the i|5.1l transit light curve model (Table[2] column 3). TTV predictions from the first (solid black, open diamonds), second (red, open 
squares), and third (blue, open circles) dynamical model in Table [3] All three models match the data well. Right: Same models as left 
plotted over longer timespan; the models differ in their predictions for future O-C variations. 



TABLE 3 

Parameter values for TTV fits. Fixed in all fits are A'h = 1.22Af Q , M.oi = lM Jup , i.oi = 90°, and O.oi = 0°. Orbital 

ELEMENTS ARE JACOBIAN ELEMENTS (THE OUTER BODY'S ORBIT REFERRED TO THE CENTER-OF-MASS OF THE STAR AND THE PLANET) 

DEFINED AT DYNAMICAL EPOCH BJD 2455200. 
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0.74 


90° 


660.7 


496.0 


-0.0092 


-0.1824 


6.66 


90° 


0° 


4.65 


±0.001696 


±0.002734 


fixed 


fixed 


±21.0 


±7.2 


±0.0105 


±0.0192 


±0.34 


fixed 


fixed 




69.721695 


238.150714 


0.90 


180° 


643.8 


304.81 


0.148 


-0.0496 


5.82 


90° 


0° 


2.62 


±0.002548 


±0.004422 


fixed 


fixed 


±50.6 


±2.31 


± 0.059 


±0.0103 


±0.98 


fixed 


fixed 




69.749706 


238.303853 


0.74 


90° 


1038.0 


841.9 


-0.0681 


-0.3567 


24.28 


60° 


130° 


0.01 


±0.000499 


±0.000672 


fixed 


fixed 


±38.5 


±21.3 


±0.0078 


±0.0148 


±0.41 


fixed 


fixed 





orbits allow for a more massive perturbing companion. 
This fit (Table [3l row 5-6), featuring a 24.3 Mj up brown 
dwarf companion with a one-thousand day orbital period 
and moderate eccentricity, is an excellent match to the 
observed TTVs and is plotted in Figure [9] In contrast 
to the coplanar fits, this model predicts deviations not 
only in the centr al transit times but in the duration of 
the t ransits (e.g. iMiralda-Escudel [2002t iNesvornv et al.1 
12012ft . due to a secular variation in the transiting planet's 
duration. However, the small transit duration variations 
predicted by this model would not be significantly de- 
tected in the current data and, depending on the impact 
parameter, may or may not be detectable in by the Ke- 
pler extended mission. Comparing the goodness of this 
fit to the two coplanar ones, we see that we can neither 
distinguish the orbital plane of the third body, nor limit 
its mass to the planetary regime. 

In all three cases, we see in the integrations that, as ex- 
pected, the "jumps" in the TTVs correspond to the com- 
panion's periapse passage. In the right panel of Figure[9l 
we plot the TTVf0 into the future. Additional transits 
in the Q7-Q12 data scheduled for future public release 

14 In plotting these extended models, we have slightly adjusted 
the linear ephemeris of the transiting planet to remain consis- 
tent with the data while keeping future O-C variations centered 
at 0. Otherwise the predicted differences between the three differ- 
ent models appear misleadingly large. 



and through the Kepler extended mission may allow us 
to distinguish among them, as well as the many other 
possible models among which we cannot distinguish cur- 
rently. We have used the Bulirsch Stoer integrator in 
Mercury ( Cha mbers! 119991 ) to confirm that all three fits 
described here are dynamically stable over 10 Myr, with 
no planet-planet scattering occuring during this interval. 
The fits do not rule out past planet-planet scattering: 
in the context of HEM, the bodies could have under- 
gone scattering in the past and subsequently stabilized as 
KOI-1474.01's orbit shrank through tidal dissipation. We 
note that the transiting planet's eccentricity undergoes 
secular variations and, in the case of the first two fits, the 
current eoi is not the maximum and thus the planet ex- 
periences enhanced tidal dissipation during other parts 
of the secular cycle. We discuss this behavior further 
in the next section, in which we consider whether KOI- 
1474.01 is a failed- or proto-hot Jupiter. We defer ex- 
haustive exploration of the parameter space of the three 
body model until more data are available, including ad- 
ditional transit times that extend the baseline to cover 
the perturber's subsequent periapse passage and comple- 
mentary constraints on the perturber's mass, period, and 
eccentricity from planned radial-velocity measurements. 
However, the possibilities illustrated here show that pin- 
ning down the perturber's mass and orbit will likely re- 
veal clues about the past mechanism of HEM and the 
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future fate of KOI-1474.01. 

6. KOI-1474.01: A PROTO- OR FAILED-HOT JUPITER? 

KOI-1474.01 is a highly eccentric, Jupiter-sized planet 
being perturbed by an unseen companion, the "smoking 
gun" that may have been responsible for KOI-1474.01's 
HEM. The transiting planet might be either a proto-hot 
Jupiter that will achieve a short period, low eccentric- 
ity orbit via tidal dissipation over its host star's lifetime 
or a failed-hot Jupiter, too far from its star to experi- 
ence significant tidal dissipation. If the planet is a failed- 
hot Jupiter, it is destined to spend the rem ainder of its 
host star's lifetime in the "period valley" dJones et al.l 
[20TjllUdrv et alJl2003tlWittenmver et al.l2010t l. between 
the region where it formed (beyond 1 AU) and the hot 
Jupiter region (P < 10 days ks 0.091 AU). 

S12 predicted the discovery of super-eccentric hot 
Jupiter progenitors among the Kepler candidates based 
on the following argument. A Jupiter kicked to a small 
periapse via one of several proposed HEM mechanisms 
will enter the proto-hot Jupiter stage. Assuming that a 
steady flux of hot Jupiters are being spawned through- 
out the Galaxy, there must exist a steady-state stream 
of highly eccentric planets on their way to becoming the 
population of hot Jupiters thus far observed. The tidally- 
decaying Jupiters follow tracks of constant angular mo- 
mentum: Pfi„ai = P(l — e 2 ) 3 / 2 , where P and e are the 
values corresponding to any time during the circulariza- 
tion process. 

To predict the number of highly-eccentric proto-hot 
Jupiters that Kepler will discover, S12 used the Exo- 
planet Orbit Database (EOD) sample of planets with 
M„shu > 0.25Mj up and Pfi na i < 10 days (|Wright et al.l 
2011, http://www.exoplanets.org). The Pfi na i cut-off 
is motivated by the excess of currently known Jupiter- 
mass planets on circular orbits with P < 10 days. 
They computed the fraction of Jupiters in the ranges 
3 < -Pfinai < 5 days and 5 < Pfinai < 10 days that are 
moderately eccentric (0.2 < e < 0.6). Next they multi- 
plied these fractions by the total number of Jupiter-sized 
(R > 8i?®) Kepler candidates in these two Pfinai ranges, 
yielding the predicted number of mo derat e ly ecc entric 
Kepler Jupiters. Finally, they use the Eul dl98l tidal 
equations to compute the relative number of highly ec- 
centric to moderately eccentric Jupiters at a given P nna i 
and predict 5-7 super eccentric Jupiters in the Kepler 
sample with e > 0.9 and P < 93 days. 

Because of the uncertainty in KOI-1474.01's eccentric- 
ity, we cannot definitively say whether it is one of the 
super-eccentr ic Ju piters predicted by S12. From our or- 
bital fits in £|5.1[ we derive a P nn ai posterior distribu- 
tion of which 42% have Pfinai < 10 days and 19% have 
Pfinai < 5 days. Therefore, the evidence only slightly fa- 
vors the interpretation that KOI-1474.01 is a failed-hot 
Jupiter with a P nn ai > 10 days. Follow-up, high-precision 
radial velocity measurements may allow us to constrain 
KOI-1474.01's eccentricity even more tightly and con- 
firm or rule out e > 0.9 and P nna i < 10 days. Further- 
more, the perturbing companion may cause secular vari- 
ations in KOI-1474.01's eccentricity ( ^5.3|) . boosting the 
tidal circularization rate during intervals of higher eccen- 
tricity; additional constraints on the perturber's identity 
may one day allow us to explore this effect. 

In Figure [lOl we display KOI-1474.01 (gray circle) in 



the context of the current sample of Jupiter-sized and 
Jupiter-mass planets. We plot the quantity (1 — e 2 ) 
vs. a to allow us to overlay tracks of constant an- 
gular momentum while visually distinguishing high vs. 
low eccentricities. An ag na i track is the path through 
phase space that a particular Jupiter follows during its 
tidal evolution; a Jupiter's current a nna i defines its an- 
gular momentum and remains constant as the Jupiter 
undergoes tidal circularization. The solid, black lines 
represent tracks of angular momentum corresponding to 
flfinai = 0.057,0.091 AU, i.e. Pfinai = 5,10 days around 
Sun-like stars. Any Jupiter along an afi na i track will stay 
on that track, reaching a = a nna i as its e — > 0. The other 
symbols represent planets with M p sini > 0.25Mj up , 
0.7M© < M* < 1.3M(7), and m easured eccentricities 
from the EOD (j Wright et alJl201lT) . The median of KOI- 
1474.01 's eccentricity posterior places the planet in the 
period valley from 0.1 < a < 1 AU, along with 
about a dozen other eccentric Jupiters. At one-sigma, 
KOI-1474.01 may be within (i.e. to the left of) the 
afinai < 0.057 AU track (i.e. will end up at a semi- 
major axis less than 0.057 AU if it fully circularizes), 
like the poster-planet of high eccentricity, HD 80606 b 
(red square). 

However, KOI-1474.01's ultimate fate is determined 
not only by P nna i but by its tidal circularization rate; 
even if the planet has Pfinai < 10 days, it will not be- 
come a hot Jupiter unless it can circularize over its 
host star's lifetime. A hot Jupiter's tidal circulariza- 
tion rate depends on a combination of orbital properties 
an d physical planetary a nd s tellar properti es. Follow- 
ing |EggletoneriI] (|1998[ ) and|Hansen| (|2010[ ) — and ne- 
glecting the effects of the planet's spin and tides raised 
on the star — a tidally-circularizing planet's eccentricity 
e-folding time is: 



a 8 (l 
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where ap is the planet's internal dissipation constant and 
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Note that the tidal circularization timescale e/e depends 
steeply on the planet's semimajor axis and eccentricity, 
but only weakly on physical stellar and planetary param- 
eter^. Therefore we might expect to see a signature of 
tidal circularization in our 1 — e 2 vs. a plot even ne- 
glecting the difference in physical properties among the 
planets plotted. 

First imagine if all the planets underwent HEM at once 
and have tidally evolved for time t. A certain curve in 
(1 — e 2 ) vs. a space, a c i rc (a, 1 — e 2 ), represents the cir- 

15 The other parameter raised to a large power is Rp° ■ Most 
objects with M > 0.25Aifj up — from Jupiters to brown dwarfs 
— have R p ^ Rj U p\ the R p ° term varies by a factor of 60 from 1 
Jupiter radius to 1.5 Jupiter radius. However, in practice wc find if 
that w e no rmalize a by (R p / Rj U p) 5 ' 4 for planets with known radii, 
Figure nol does not change significantly. The circularization rate's 
strong dependence on a dominates, because a undergoes large frac- 
tional changes throughout the hot Jupiter region, with a change in 
semimajor axis of 0.02 AU corresponding to an order of magnitude 
change in the tidal circularization timescale. 
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Fig. 10. — Distribution (1 — e 2 ) vs. semimajor axis for known exoplanets from the EOD with 0.7Mq < Ai* < 1.3Mq, measured 
eccentricities, M p sini > 0.25Mj up , and with apoapses beyond 0.9 AU (blue, downward triangles) or within 0.9 AU (red, upward triangles) 
(Wright et al. 2011, http://www.exoplanets.org). The gray circle marks KOI-1474.01, with the asymmetric gray error bars representing 
the 1-sigma (dark gray), 2-sigma (light gray) confidence interval of KOI-1474.01's eccentricity. HD 80606 b is denoted with a red square 
symbol. The solid black lines are tracks of constant angular momentum corresponding to afj na i = 0.057, 0.091AU; each indicates a track 
that a single Jupiter follows through phase space as it undergoes tidal circularization and maintains a constant angular momentum. As 
it fully circularizes (e — > 0), a Jupiter ends up at the top of the track at 1 — e 2 = 1. The purple, dashed lines represent constant tidal 
circularization rates corresponding to a c j rc = 0.057, 0.091AU (Eauation ll3jl , A group of Jupiters that lie along a particular purple a c i rc 
line is undergoing tidal circularization at the same rate. 

cularization time (Equation [TTj) equal to t. We would 
expect this curve to envelope the still-eccentric Jupiter 
population, because all planets to the left of the curve 
(i.e. with 1 — e 2 less than the curve for a given semi- 
major axis) would have already undergone an e-folding's 
worth of circularization. The semi-major axis a = a c ; rc 
would be the edge of the circular population we call "hot 
Jupiters," planets for which t was a sufficient amount 
of time to circularize. In reality, proto-hot Jupiters 
are being continuously spawned as new stars are born 
and as Jupiters undergo HEM. However, because of the 
steep a tidal dependence - with the tidal circularization 
timescale changing by an order of magnitude roughly ev- 
ery 0.02 AU in the hot Jupiter region — we still expect 
to see an a c ; rc boundary, corresponding to a circulariza- 
tion time equal to a typical stellar lifetime. To the left 
of this this a c ; rc boundary would be only true proto-hot 
Jupiters, caught in the act of tidal circularization. With 
a detailed accounting for observational bias and the rela- 
tively weak effects of the planets' different physical prop- 
erties, one could predict the rela tive number of proto-hot 
Jupiters on each a c i rc curve (e.g lHans"enll2010[ ). 

Solving Equation (fTT|) for (1 — e 2 ), we can combine all 
the constants — including the timescale e/e — into a c j rc 
and rewrite: 



(l- e 2) /e -3/13=(^ ) 



16/13 



(13) 



where a c i rc represents the distance within which circu- 
lar hot Jupiters have arrived via tidal dissipation. For 
small eccentricities, the factor of f e is negligible. For 



large eccentricities, we can solve Equation (|13p numer- 
ically for (1 — e 2 ). We plot a C j rc curves - along which 
all Jupiters have a similar tidal circularization rate - in 
Figure [10] as purple dashed lines. We emphasize that al- 
though the black afi na i lines and purple, dashed a C j rc lines 
Figure[in]are close together, their physical interpretation 
is different: the quantity a c i rc represents a proxy for the 
tidal circularization rate, whereas afi na i is a track that 
an individual Jupiter follows as it undergoes tidal circu- 
larization obeying conservation o f angular mom e ntum . 
If the tidal evolution according to lEggleton et al.l ()1998[ ) 
that yielded Equation [13] is a good approximation, then 
icirc may be the best quantity to consider for the cut-off 
between proto- and failed-hot Jupiter. 

Since we see a pile-up of circular hot Jupiters and 
no Jupiters with 1 — e 2 < 0.9 to the left of the pur- 
ple dashed line a c j rc < 0.057 AU (P = 5 days around a 
Sun- like star), this may represent the timescale at which 
circularization happens over a fraction of a stellar life- 
time. Under this interpretation, HD 80606 b's identity 
as a proto-hot Jupiter is not certain: it lies between 
flcirc = 0.057 AU and a C j lc = 0.091 AU, along with sev- 
eral other eccentric Jupiters that have yet to circular- 
ize. Using the internal dissipation c onstant op = 3 .4 x 
10~ 7 (5.9 x 10" 54 )g" 1 cm- 2 derived bv lHansenl (pOlOl ). the 
cut-off is even stricter: a Jupiter-like planet around a 
Sun-like star would only undergo an e-folding's worth of 
circularization over 10 Gyr if i t had a c j rc < 0.034 AU. 
However, we note that iHansenl (|2010ft derived the tidal 
dissipation constant under the assumption that proto- 
hot Jupiters, upon beginning their tidal circularization, 
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have eccentricities drawn from a normal distribution with 
a mean e = 0.2 and standard deviation of 0.25. If the 
starting eccentricities are larger — as assumed by S12 for 
proto-hot Jupiters — a larger dissipation constant may 
be necessary to match the observed hot Jupiter sample. 
In order for a 10 Gyr e-folding time to correspond to 
Qcirc = 0.057 AU, the dissipation constant would need to 
be larger by a factor of 60. 

The two-sigma upper limit on KOI-1474's eccentricity 
places the planet within a c i rc < 0.057 AU, but the two- 
sigma lower limit places it well beyond this boundary. 
The host star's age is currently poorly constrained 
f fc|3.3[) . and we do not know how recently the planet 
underwent HEM. However, if the assumptions behind 
the discussion above are correct, the steep dependence 
of the tidal circularization rate on a and e means that 
most Jupiters within a c i rc < 0.057 AU would have cir- 
cularization timescales <C and most Jupiters beyond 
icirc > 0.057 would have circularization timescales 3> t*. 
Thus the planet's fate is not sensitively dependent on ei- 
ther the star's age or when the planet underwent HEM; 
the more important quantity to pinpoint is e. 

Finally, we note that the expected number of proto-hot 
Jupiters depends on the timescale for the S12 assump- 
tion of steady production. Consider the following two 
possibilities for the dominant HEM mechanism: 

• HEM typically occurs on a short timescale com- 
pared to the stellar lifetime (for example, imme- 
diately as the gas disk has dissipated). Since we 
cannot detect planets via the transit or radial- 
velocity method around very young stars due to 
their enhanced activity, we would miss most proto- 
hot Jupiters, except for those in the small sliver 
of parameter space for which tidal circularization 
timescale is of order one stellar lifetime. 

• HEM typically occurs on a timescale comparable to 
the stellar lifetime. In this case, we would expect 
to see proto-hot Jupiters at every a C i rc , with the 
relative number of eccentric Jupiters (accounting 
for observational biases) set by the tidal circular- 
ization timescale corresponding to that a c i rc . 

The timescale of HEM depends on which HEM mech- 
anism is at play and on the typical initial architectures 
of planetary systems (e.g. for the planet-planet scatter- 
ing mechanism, how tightly packed the initial configura- 
tion is). Therefore, the discovery of definitive proto-hot 
Jupiters would not only reveal that HEM occurs but also 
constrain the details of the dominant HEM mechanism. 
If the highly eccentric planets we find are clustered at a 
single a c i rc — which would correspond to a tidal circu- 
larization timescale of order the stellar lifetime — then 
we would conclude that HEM usually occurs early in a 
planetary system's history. But if highly eccentric plan- 
ets are found at a range of a c ; rc — including a C i rc within 
(i.e. to the left of) which most planets have circularized 
- then we would conclude that HEM typically occurs 
throughout a planetary system's history. 

7. DISCUSSION AND FUTURE DIRECTIONS 

We have identified KOI-1474.01 as a highly eccentric, 
Jupiter-sized planet using a combination of a detailed 



analysis of the light curve sh a pe an d the statistical val- 
idation procedure of iMortonl (|2012l ). This makes KOI- 
1474.01 the second planet or planet candidate with an ec- 
centricity measured solely via the duration aspect of the 
"photoeccentric effect," joining KOI-686.01 whose eccen- 
tricity we measured in Paper I. We measured one com- 
ponent of the angle between the stellar spin axis and the 
planet's orbit, finding that the degree of misalignment is 
not currently well-constrained. Based on the variations 
in KOI-1474.01's transit times, we explored the identity 
of a perturbing companion; we found the TTVs to be 
consistent with perturbations from a massive, eccentric 
outer companion but could not uniquely constrain the 
perturber's mass, period, eccentricity, and mutual incli- 
nation with the currently available data. However, the 
main reason the perturber's parameters are poorly con- 
strained is that we have only witnessed perturber peri- 
apse passage; we are likely to witness another periapse 
passage over the timespan of the Kepler mission, poten- 
tially allowing us to distinguish between possible per- 
turbers, including a coplanar giant planet vs. a brown 
dwarf with a large mutual inclination. 

Because of the uncertainty in KOI-1474.01's measured 
orbital eccentricity and possible secular variations in that 
eccentricity due to the perturbing companion, it is not 
yet clear whether KOI-1474.01 is a proto-hot Jupiter — 
with a periapse close enough to its star that the planet 
will undergo full tidal circularization over the star's life- 
time — or a failed-hot Jupiter, just outside the reach 
of fast tidal circularization. However, either way, the 
planet's discovery adds to the growing evidence that 
HEM mechanisms play a major role in shaping the ar- 
chitecture of planetary systems. Th e broad eccentric- 
ity d istribution of extrasolar planets (jJuric fc Tremainel 
2008), the sculpting of debris disk s by planets 
on in c lined and eccen t ric o rb its (e.g IMouillet et all 
19971 jThommes et all 119991: lAugereau et all | 2001 



Quillenl 120061 jLevison et all 120081: Ichiang et all 12009: 



Dawson et ail 120111: iDawson fc Murrav-Clavll2012D. th e 
population of free-floating planets ( Sumi et al.l 120 111 ), 
and the large mutual in clinations measured in the Up- 
silon Andromeda system ([McArthur et alll2010f) all point 
to a dynamically violent youth for planetary systems. 
But the strongest evidence for HEM comes from hot 
Jupiters themselves — their existence an d, in many 
cases , misaligned or retrograde orbits (e.g I Winn et all 
I2009H Johnson et alJ[20TTl : lTriaudll2011l ). 

As a proto- or failed-hot Jupiter, KOI-1474.01 plays 
the crucial role of linking hot Jupiters, which are intrin- 
sically rare, to other planetary systems. Even though 
they ma ke up only a small percentage of the planet pop- 
ulation (IHoward et al.ll20Tol: IHoward et alll2012t lYoudinl 
[Mill: IMavor et all 120111: I Wright et all \2QV$ )we focus 
attention on hot Jupiters because, like meteorites dis- 
covered in Antarctica, they are known to come from 
somewhere else, bringing with them vital informa- 
tion about the past. In contrast, we do not know 
whether planets at greater orbital distances or of smaller 
sizes underwent migration, or if they formed in situ 
(e.g. IVeras et al.l 120091: fHansen fc Murray! I2012D . More- 
over, the HEM mechanisms for produ cing hot Jupiters 
- in cluding planet-planet scatte ring dNagasawa fc Id"al 
20111). the Kozai mechanism dWu fc Murray! 120031: 



Fabrvckv fc Tremainel 120071: iNaoz et all 1201 ll ) . dvnam- 
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ical r elaxation (lJuric fe Tremaind (2008), and secular 
chaos (Wu & Lithwick 2011) — make specific predictions 
for the inclination distributions of hot Jupiters, which 
can be probed via the Rossiter-McLaughlin effect. The 
existence of proto- and failed-hot Jupiters will allow us 
to argue that the mechanisms for producing hot Jupiters 
are, more generally, the mechanisms that sculpt many 
types of planetary systems, particularly those with giant 
planets within 1 AU. 

The KOI- 1474 system — an inner proto- or failed- 
hot Jupiter with a massive, long-period companion — 
may be the prototype of systems of hot Jupiters 
with dista nt, massive, outer companions , including as 
HAT -P-13 (jBakos et alJl2009h . HAT- P-17 (IHoward et all 
2012t a hot Saturn ), and Qatar-2 (jBrvan et al.l 120121) . 
Bryan et al.l ()2012[ ) present a compilation of the eight 
other hot Jupiters with kno wn outer companions. 
HD 163607 (|Giguere et alJl20"H resembles KOI-1474.01 
in that it harbors both an eccentric inner planet (e 
= 0.73, P = 75.29 days) and an outer companion 
(in this case, a massive outer planet); however, inner 
planet HD 163607 b is very likely a failed-hot Jupiter, 
as it has Pfi na i = 24 days. The expanding baseline 
for radial-velocity measurements may reveal additional, 
long-period outer companions of other hot Jupiters, 
proto -hot Jupiters, and failed-hot Jupiters (jWright et al.l 
2009). These additional companions may have been the 
culprits responsible f or the HEM o f their inner brethren. 
Moreover, although Steff en et al.l (|2012f ) examined the 
transit timing variations of Kepler hot Jupiters and 
found no evidence for nearby massive planets, the ex- 
tended Kepler Mission will allow for the detection of 
distant companions, should they exist, through TTVs. 

Through radial-velocity follow up with Keck/HIRES 
we will measure the mass of KOI-1474.01, tighten the 
measurement of its high eccentricity, place additional 
constraints on the outer companion, and potentially dis- 
cover additional bodies in the system. Assuming a 
Jupiter-like composition to estimate a mass for KOI- 
1474.01 of M p w M Jup , host star KOI-1474 would have 
an radial velocity semiamplitude of ~ 70 m s _1 , feasi- 
ble for detection using Keck/HIRES. We will then com- 
bine the RV-measured eccentricity with the transit light 
curves to more tightly constrain the stellar parameters, 
yielding a better constraint on the planet's line-of-sight 
spin-orbit angle \i — i s \ } which is currently ambiguous due 
to uncertainty in the stellar radius. It may even be pos- 
sible to detect the Rossiter-McLaughli n effect, whic h has 
a maximum amplitude of ~ 50 m/s ()Winnl 1201 Ol eqn. 
40). Although RV measurement s of such a fai n t star 
(Kp = 13.005) pose a challenge, Uohnson et al.l (|2012[ ) 
have demonstrated the feasibility of following up faint 
Kepler targets with their measurements of KOI-254, a 
much fainter, redder star (Kp = 15.979). 

KOI-1474.01 contributes to the growing sample of 



proto- and failed-hot Jupiters. From an estimate of the 
unbiased number of proto-hot Jupiters, we can determine 
whether HEM accounts for all the hot Jupiters observed, 
or whether another mechanism, such as smooth disk mi- 
gration, must deliver som e fraction of hot Jupiters. (See 
iMorton fc Johnsortll2011al for the statistical methodology 
necessary for such a measurement.) Transiting failed- 
hot Jupiters orbiting cool stars will be v aluable tar- 
gets f or testing the obliquity hypothesis of I Winn et al.l 
(2010) that hot Jupiters realign cool stars: we would ex- 
pect failed-hot Jupiters - which have long tidal friction 
timescales — to be misaligned around both hot and cool 
stars. 

Designed to search for Earth twins in the habitable 
zones of Sun-like stars, Kepler is revealing a wealth of 
information about the origin of the most unhabitable 
planets of all: hot Jupiters. Kepler's precise photom- 
etry, combined with a loose prior on the stellar density, 
allow us to measure the eccentricities of transiting plan- 
ets from light curves alone and to search for the highly 
eccentric proto- and failed-hot Jupiters we would expect 
from HEM but not from smooth disk migration (S12). 
If our basic understanding of HEM and tidal circulariza- 
tion is correct, KOI-1474.01 is the first of a collection of 
highly eccentric planets that will be discovered by Kepler 
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